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ABSTRACT 

The exact profile of a gamma-ray burst (GRB) afterglow image on the plane of 
the sky can provide important constraints on the underlying physics. In particular, it 
can test whether the magnetic field in the emitting shocked external medium decreases 
significantly with the distance behind the shock front, or remains roughly constant. 
Moreover, it enables more accurate measurements of the afterglow image size and the 
expected scintillation properties. In this work analytic expressions are derived for the 
afterglow image in power law segments (PLSs) of the afterglow synchrotron spectrum 
in which the emission originates from a very thin layer just behind the shock, while 
simple semi-analytic expressions are derived for the remaining PLSs in which the 
emission arises from the bulk of the shocked fluid. In all cases the expressions are for 
a general power law external density profile, and are convenient to use in afterglow 
studies. 
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1 INTRODUCTION 

The direct measurement of the size of the gamma-ray 
burst (GRB) afterglow image in the rad io, both through 



the quenching of diffractive scintillation (Fra il et alj 



and more directly us i ng th e VLB A (Taylor et al.l 



1997) 



2004 



less certain and robust. For this r eason, most ca l cula- 
tions of the afterglow image so far JWaxman 1997 : Sari 
19981: iPanaitescu fc Meszarod [l998l : iGranot. Piran fc Sari 



\ l2005l : [ Pihlstrom et al.l l2007h . provided good support for 
the basic dynamical pictur e of standa rd afterglow the- 
ory (Waxman. Kulkarni fc F rail 1998; Ore n. Nakar fc PiranI 



I2OO4I : IGranot. Ramirez-Rui z fc Locb 2005'). The surface 



brightness profile within the afterglow image can provide 
additional constraints on the afterglow physics. It can poten- 
tially be tested directly for a particularly nearby and reason- 
ably bright GRB afterglo w, or even for more distant events 
in the case of microlensing (iGarnavich. Loeb fc Stanekl2000l : 
iGaudi. Granot fc Loeb|[200ll . e.g.). In particular, the exact 
appearance of the afterglow image on the plane of the sky 
can be very useful in improving the accuracy of the after- 
glow image size measurements, both when the image is di- 
rectly but only marginally resolved and its size is determined 
through fits to the visibility data, and through the quench- 
ing of scintillation. It can also improve the estimates for the 
expected amplitude of scintillation, and thus help in the af- 
terglow modeling. 

The dynamics of GRB afterglows before the jet break 
time, tjst , are well described by the iBlandford fc McKed 
j 19761 . hereafter BM76) self-similar solution, while at later 
times the exact dynamics of the GRB jet are much 



1999al lbl: IGranot fc Loeb 20011 ) have been for the (quasi 
) spherical stage, corre sponding to t < tjct (see, however, 
lloka fc Nakamural 120011 ). We shall also address the (quasi- 
) spherical stage, for the same reasons. In previous works 
the expressions for the afterglow image are either not very 
accurate due to some simplifying assumption, or not fully 
analytic, which makes them inconvenient for others to use. 
Therefore, in § [2] we derive analytic (§§ 12.11 12.2|l or semi- 
analytic (§ 12. 3p expressions for the surface brightness dis- 
tribution within the afterglow image, for a general power 
law external density, and for all of the power law segments 
(PLSs) of the afterglow synchro tron spectrum that are de- 
scribed in IGranot fc Saril (|2002| . hereafter GS02), using the 
notations of GS02. The magnetic field is assumed to be tan- 
gled on small scales with an isotropic distribution in the co- 
moving frame of the emitting shocked fiuid, and hold a con- 
stant fraction of the internal energy everywhere. The main 
results are shown in Figures [T] and [51 

2 DERIVING ANALYTIC EXPRESSIONS FOR 
THE AFTERGLOW IMAGE 

Consider a spherical relativistic blast wave. A photon that 
is emitted from the shock front at a lab frame time t when 
the shock radius is R and from an angle 6 relative to the line 
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of sight to the central source (that is located at the origin 
reaches the observer at an observed time 

R 



tohs — i 



■ COS 6 



(1) 



where tohs = corresponds to a photon emitted at the origin 
(i.e. the central source) ai t = [i.e. the time when the 
outflow was launched, R{t — 0) ~ 0]. For convenience we 
normalize the shock radius by its maximal value along the 
equal arrival time surface of photons to the observer (along 
the line of sight, at = 0), y = R/Ri- For a power law 
external density profile, poxt = Ar~^ with fc < 4, the Lorentz 
factor of the shock front during the relativistic phase (F 3> 
1) scales as F oc rC'-^^'^ (BM76), i.e. F = Tiy^'"'^^^^. This 
implies 



c 



f 

Jo 



dR 



2cF2 2(4-fc)Ffc 



, 4 



(2) 

and that Ri = 2(4 - k)Tictahs = 4(4 - fc)7fctobs, where 
7; = F;/\/2 is the Lorentz factor of the fluid just behind the 
shock a,t R — Ri. Since we are interested in the relativistic 
regime (F 3> 1), this implies that all the relevant emission 
is from small angles {9 <^ 1) so we can approximate cos 9 ~ 
1 — 9'^ /2, and equations [T] and [2] imply 

,,4-fc 



1 



y 



The distance from the center of the image is given by 

Ri 



Ri 



Rsm9 : 



V: 



y-y 



and its maximal value is 

R± ,max 



■(5-fc)/2(4-fe) 



(3) 



(4) 



(5) 



and is obtained at y, = R*/Ri = (5 — fc) and 
9, = Fr' = F-i(5 - fc)-(3-'=)/2{4-fc)^ ^ljgj.g ^ p(^^)^ 

Therefore, the normalized distance from the center of the 
circularly symmetric afterglow image is given by 



R± 



(5-fc) 



(5-fc)/2(4-fc) 



Ri 



■yy- 



(6) 



2.1 Self- Absorbed Spectral Power Law Segments 

Below the self-absorption frequency !/sa (PLSs A, B, C in 
GS02) the specific intensity (i.e. surface brightness) 1^ is 
equal to the source function [Sv = jv/otv) at the front end 
of the equal arrival time surface, (5 — A:)~^'''*~*' ^2/^1 
{R, ^ R^ Ri), where y = (5 - k)-^'^^'''^ (R = R,) and 
y = 1 {R = Ri) correspond to 2: = 1 and x = 0, respectively: 



(7) 



where I^, cc (i^')*, kTcs{i^) is the effective temperature of the 
electrons that radiate in the local rest frame of the emitting 
fiuid just behind the shock at the observed frequency v, and 



27 



2^/2(4 -fe)F,y(5-'='''2 



1 



(8) 



-726/2 (7 - 2fc)2/4-'= + 1 

is the Doppler factor. Primed and un-primed quanti- 
ties are measured in the comoving (emitting fluid) and 
lab (or observer) frames, respectively. In PLS A {i/m < 



V < 

^syn ("T* 

,.i/2r- 



i^sa), kTcB{i') ~ 7emeC for 7e that satisfles v « 
; eB'7^/(27rmec), so that kT^a oc {v/B'^^^ oc 



-i/2„-i/4 



Pcxt oc U 

and 6 = 5/2. Thus 



l/2j^3/4 



[since {B') oc ej^t oc pcxtF 



5/2 k/4 

V ' y ' 



1 + 



2(4-fe)y4 



-1/2 



(9) 



h{x = 1) 
h{x = 0) 



'0.8165 (fc = 0) 



_(5 _ fc)4r4^ 



1^0.6204 (fc = 2) . 

In PLS B (i^ < min[i^m, !^sa, I'ac]), kTcH ~ "imrrieC? oc F oc 
y(k-z)/2^ and therefore 6 = 2 and 



2{4-k)y 



(10) 



hjx = 1) 
L{x = 0) 



fe)' 



'2.2291 (fc = 0) 
1.1547 (fe = 2) 



In PLS C (!^ac < I' < Vsa) the emission is from elec- 
trons that have suffered signiflcant cooling. Locally, the 
electron distribution in this region is approximately mono- 
energetic, and the electron Lorentz factor scales as 7e oc 
1/[{B')^1'] with the distance /' behind the shock front (at 
the comoving time of emi ssion) in the comoving frame 
l|Granot. Piran fc Saril2000l) . Most of the photons that reach 
the observer near an observed frequency f originate near 

n') which is given by t,^[Ii {f, h')] = 1, where n' = cos 9' , 
and 9' is the angle between the direction of the photon and 
the shock normal (i.e. the radial direction) in the comov- 
ing frame. Therefore, kTeeiv) « ■ye[li{i', tJ-')]meC^ . The path 
length of a photon until it overtakes the shock front (which is 
receding from the shocked fluid at a velocity of /J'^i^c ~ c/3), 
s' , is related to its initial distance from the shock front, I' , 
by 

s' _ 1 _ 3(1 -h7^6|2 
T ^ 



/ 



^i' - 2(1 - 27292) 



(7-2fc)y* 



+ 1 



(5 - k)y* 



- 1 
(11 



- 7 



-Pfi 1 + 72612 

The location of I'l is where the optically thin and optically 
thick fluxes are equal, 

P', / „' \i/3 



n's'(Z'i)- 



1^' ,max 



47r 



7e(<i)meC 



(12) 



oc B'jj oc R ^'^^7e, n' oc 



where P^,^^^^ oc B' oc R'^^^, 
Fpcxt oc _R-{3+fe)/2^ and -yeH'i) oc l/[{B'fl[]. In order to use 
eq. (O we evaluate eq. (|12p at u' = so that oc 

Z/5/8/-3/8(^')-3/8(5,)-3/2 ^y^^(^) ^ le[l'l{y)] OC 



6 = 11/8 and 



-1/2 



5/8 j3/8y(3-3fc)/16 



implying 



11/8 (5fe-18)/8 

(X u ' y^ " 



-13/8 



2(4 - fc)y4-fc 
(7-2%^-'= 



+ 1 



(5 - k)y^ 



3/8 



(13) 



Note that the surface brightness diverges at the outer edge 
of the image (x = 1), as oc (1 — x)~''^^^^ for 1 — x ^ 1. 
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2.2 Fast Cooling Spectral Power Law Segments 

For PLSs F (max[!/c, J^sa] < f < Vm) and (v > 
tlie emission comes from electrons that cool 



on a time-scale much smaller than the dynamical time, and 
therefore it originates from a very thin layer just behind the 
shock so we can use the values of 7, n' and B' just behind 
the shock. We have 



.^'jt' ~/«c[7e(^')]j^ , (14) 

where / is given by eq. (HIlEI, I'^iye) ~ 2nml/[aT{B'f^e] 
is the electron cooling length, and v' = fsyn[7e(j'')] ~ 
eB'-fl{v')/{2-Kmec). Let L'^, oc R^iv'f be the spectral lu- 
minosity (the total emitted energy of the whole shell [i.e. 
thin emitting layer of shock fluid] per unit time and fre- 
quency, assuming a spherical emitting shell), and P'^i be 
the energy emitted per unit time, volume, and frequency 
(where both are measured in the comoving frame). For 
PLS F, a = (5 - 2fc)/4 and 6 = -1/2, while for PLS 
14 - 9p -I- 2k (p - 2)]/4 and b = -p/2 (see Ta- 
G"ranotll2005h . For isotropic emission in the co- 
moving frame we have j'^, — P^;/47r and therefore L'^^, ^ 
AnR'^U'y,{u')]P;^, oc 7?^/^[7e(i'')lj^' « where both 

P^/ and jl, axe evaluated inside the thin layer of width 
l'c["/e{j^')] behind the shock front in which the electrons 
whose synchrotron frequency is v' have not yet cooled signif- 
icantly (and L^, is evaluated using the volume of this layer). 
Therefore we have /,y oc fR^'^iuf and oc fR^-^u'' so 
that 



H, a = 
ble 1 of 



|(5-fc)t/4 



1| [(7- 2fc)y4-fc + 1]^ 



(15) 

In order to express as a function of x rather than y we 
use eq. ((Sjl to obtain y(x). It is important to note that 
y{x) it is double valued, where the two values correspond 
to the front (j/+ J? y*) and back (j/_ ^ ?/,) of the equal 
arrival time surface (EATS) of photons to the observer. 
Here y-{x) = R-{x)/Ri corresponds to a photon emit- 
ted from the back of the ETAS at a relatively small ra- 
dius, R-, and a large emission angle, 9 > l/r(i?_), which 
corresponds to an angle 8sh > 90° relative to the radial 
direction in the rest frame of the shock front, and there- 
fore initially lags behind the shock front. Eventually, at 
some larger radius R+, it catches- up with the shock front 
and starts moving ahead of it. From R+ onwards its tra- 
jectory coincides with that of a photon emitted from the 
front of the ETAS [y+{x) = R+{x) /Ri] at the exact place 
and time when it crossed the shock front, but at an an- 
gle 9 < l/r(_R+) that corresponds to 5sh > 90° so that 
it never l ags behind the shock front (For more details see 
Fig. 1 of iGranot. Cohen-Tanugi fc do Couto e Silval \200i . 
and the related discussion therein). The two values, y-{x) 
and y+{x), coincide at j/, which corresponds to a; = 1, i.e. 
at the outer edge of the image, where the surf ace bright- 
ness diverges as 7,^ oc (1 - x)-'^/'^ for 1 - a; < 1 l|Sarilll998l : 
iGranot fc Loebll200ll ). For fc = 3 the shock Lorentz factor 



^ With the exception that here the absolute value of the denom- 
inator should be used, as is becomes negative for y < (at the 
back of the equal arrival time surface) since the photons (initially) 
move away from the shock in this case. 



does not change with radius, the equal arrival time surface 
becomes an ellipsoid, and there is a particularly simple so- 
lution: J/-I- = i (1 ± vl — 2?) . For the physically interesting 
case of fc = 2, which corresponds to the stellar wind of a 
massive star progenitor, we also obtain an explicit analytic 
solution; 



V±(x) — — ;= COS 



TV =p arctan \/ a;"^ — 1^ 



(for fc = 2) . 

(16) 

The total value of Iv{x) is obtained by summing these two 
contributions (in §[27l]only the value corresponding to y+ ^ 
y, should be used, since the back of the EATS is obscured). 

For sufficiently large values fc [fc > 32/9 ~ 3.556 for PLS 
F, and fc > 4(9-p)/(10-p) for PLS H] the surface brightness 
diverges at the center of the image {x <C 1) due to contribu- 
tions from small radu (y < 1), as h oc a:(32-9*)/2 p^g 
F and as h oc j:lMQ-p)-(io-p)k]/2 p^g jj Physically the 
divergence is avoided due to the break down of some under- 
lying assumption, e.g. for fc > 3 the shock accelerates and 
was initially non-relativistic at some radius -Rnr correspond- 
ing to j/NR = Rnr/Ri which introduces a cutoff at a;NR ~ 
(yNR/Cfe)i/2 « 1, where = (4 - fc)(5 - fc)-(5-'=)/(4-fe). 

The expression w e obt a in for Ii,{x) is slightly different 
from that obtained bv lSaril (|l998l ) for fc = 0. The difference 
arises since he did not account for the fact that the fluid 
just behind the shock moves at a different velocity than the 
shock front itself. We can reproduce his results by replac- 
ing / in the expression for 7„ with 1/| cosSahj where 9sh is 
the angle between the direction to the observer (i.e. that of 
the emitted photons that reach the observer) and the shock 
normal measured in the rest frame of the shock front (which 
moves at a velocity of /3shC ~ c/3 relative to the comoving 
rest frame of the shocked fluid), i.e. 



/ 



1 



cos f sh 



\i-r^9^\ 



ll - 272 



1(5 - k)y'^''= - 1| 
(17) 

2.3 Spectral Power Law Segments Originating 
from the Bulk of the Shocked Fluid 

In the PLSs that have been treated so far the emission arises 
from a very thin layer of shocked fluid just behind the shock 
front. This has enabled the use of the values of the hydrody- 
namic quantities just behind the shock, and simplifled the 
derivation of the surface brightness distribution within the 
afterglow image, Iv{x), resulting in fully analytic expressions 
for it. The emission in such PLSs does not depend on the 
hydrodynamic profile of the shocked fluid downstream of the 
shock transition, and responds relatively quickly to changes 
in the external density (although the fact that the contribu- 
tions to any given observed time are from a wide range of 
radii s till causes significant s moothing of the observed light 
curve; iNakar fc Granotll2007i ). 

Now we turn to calculate Iv{x) for PLSs D, G, and E, in 
which the emission originates from the bulk of the shocked 
fluid. In these cases one must specify the values of the hydro- 
dynamic quantities everywhere within the region of shocked 
fluid. For this purpose we use the spherical adiabatic self- 
similar BM76 solution. In this solution the hydrodynamic 
variables depend on the self-similar variable x, and (GS02) 



Ri 



Ri 



-1/2 



y - xy 



X 



y - Ckx^ 



(18) 
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where Ck = (4 - fc)(5 - as well as 

7 



(19) 
(20) 
(21) 



and 



XV 



{4-k)y ' " 2(4 - k)xy*-'' 
27 23/2(4-fe)ra'/'2/''-*^^' 



1 + 726*2 



(7-2fc)xj/4-fc + l 



(22) 



Assuming isotropic emission in the comoving frame, j'^i = 
Pli /47r, and using the fact that lu ~ j jvds where ju — S^jli 
and in our case ds « dR = Ridy, one obtains 



I Ax) 



2{A - k)'-^ RiFf 



2{A~kfRirf 



dy 



x{y,x)y^ ''Pj. [y,x{y,x)] 

[{7^2k)x{y,x)y*-^ + lf 



^ [y~Cux^)y''Pl,{y,x) 
^[(7-2fc)(y-C,x2)+y]2 • 



Now we need to derive Pli [y, x{y, x)] and therefore 
Pl,{y,x) for PLSs D, G, and E. For PLSs D (i^,, < 



V < Vm < Vc) and G (max[i/m, !/sa] < v < 



P', 



n'P;.^.^Av'/^Lf where Pi, 



,max,e ~ GTrUeC^ B' / {ie) 

and = t'syn(7m) « eB'7^/(27rmec) where 7™ = 

geee7(n'meC^) and g = [p - 2)/{p - 1) for p > 2. Thus, 
for PLSs D where b — 1/3 and G where b — (1 — p)/2 we 
obtain 

P' ,/^[13fc-47+K13+fc)]/[6(4-fc)] [K4-fc)-6-fe]/2 



[(7^2k)x.y'-'' + lf 



(24) 



hccv dy[x 



xy 



[7fc-23+f)(13+fc)l/[6(4-fc)l 
[6(4-fc)+4-3fc]/2 



[(7-2fc)xy*-* 



+ 1 



(25) 



where the limits of integration over y are the appropriate 
roots of the equation xiy^x) — 1, i.e. y — y^~^ = CkX^ (see 
subsection I2.2|l . 

For PLS E (!/sa < V < Vc < Vm), the emission is domi- 
nated by regions where all of the electrons have cooled sig- 
nificantly and their energy distribution is practically a delta 
function, A''(7e) ~ n' 5{')^^ - 7max(x>y))i where 7niax is given 
by eq. A12 of GS02. Using equations A8 and A9 from that 
paper, one obtains that for a constant observed time, 

^(22-5fe)/[6(4-fc)ly(l + fc)/2 
7max(x,2y) « (^{19-2fc)/[3(4-fc)] _ ;^-) ' (^6) 

Therefore, P'^, « n'^'^,n.ax,e(^^7^'max)'''^ where u[^,^ = 
J^syn(7max) « cB' 'y'^.^y, / {2imieC) . Altogether, using eq. ((23l) 
we obtain 



IAx)(xu'^'' J 



dy { y 



(4-5fc)/3 



r^(19-2'=)/[3(4-fc)] 



-, 2/3 



^{18-5fc)/[2(4-fc)] 

X [l + (7-2fc)xy'-'=]"'^'| . (27) 
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Figure 1. The afterglow images: the surface brightness normal- 
ized by its average value as a function of the normalized distance 
from the center of the image, for all of the different power law 
segments (PLSs) of the afterglow synchrotron spectrum that are 
described in GS02, and for three value of the power law index k 
of the external density, where poxt oc r~*. For PLS E, the nor- 
malization of the dashed line for A; = 2 is arbitrary, since the 
total flux diverges in that case under our assumptions (see text 
for details). 

It can be showi£| that for k > 37/26 ^ 1.423, I,y{x < 1) oc 
^2{37-26fc)/9 that it diverges at the center of the image, 
while the flux in this regime is dominated by the contribu- 
tion from x ^ 1 and thus scales as Fi, oc I,j{x)xdx oc 
1/(23 - 13fc) for k < 23/13 and diverges for k ^ 23/13 ^ 
1.769. If a lower limit for the range of integration, a;niin, is 
introduced then Fi, oc x'^f^ i3fe)/9^ g^^j-gg^ ^ divergence 
of the total flux is unphysical, and does not really occur. In- 
stead, the underlying assumptions for this PLS break down 
in this regime (fc ^ 23/13), resulting in an introduction of 

1/2 

such an Xmin ~ !/mini ^^'^ P^S E no longer exists in the 
same form. A detailed treatment of the interesting behavior 
in that case is saved for a separate work. 

Similarly, in PLS D the surface brightness diverges at 
the center of the image for k > 61/26 « 2.346 as Iv{x -C 



For a; = we have = and x 



and this still ap- 



proximately holds in the region of interest here (x ~ y ^ 1, 



where xV 



1 — Ck^ /y becomes significantly different than 



1 only for y ~ y^i^ CkX^). Thus 7^ oc / di; j/2(i4-i3fe)/9 be- 
comes dominated by the lower limit of integration, j/min ~ C^x^ , 
for k > 37/26. In this case I, oc y^^-'^'"'^^'' oc x^(37-26k) /9 ^ 
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2.5h PLSG, b = (1-p)/2, k = 

2 
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Figure 2. This figure demonstrates tlie dependence of tiie af- 
terglow images on tlie power-law index, p, of the electron energy 
distribution, in the two PLSs (G and H) where it has some effect. 
The afterglow image becomes somewhat more limb-brightened as 
the value of p increases, due to the corresponding decrease in the 
value of the spectral index b. 



1) oc 2;2{6i-26fc)/9^ ^-^e flux diverges for k > 35/13 ^ 
2.692. In PLS G the surface brightness diverges at the center 
of the image foi0 fc > (32 - 4p)/(ll - p) as I^{x < 1) oc 
^i4-5fc-(4-fc){p-i)/2^ and the flux diverges for fc > (36 - 
4p)/(ll —p). Again, the flux cannot diverge in practice, and 
this is an indication that the model assumption break down 
in those regimes. 



nificant only within a thin layer (of width A <^ R/'y'^) just 
behind the shock front, and negligible furth er downstream 
from the shock (see, e.g. iRossi fc Reesll2003l ). then this will 
affect the appearance of the afterglow image. In particular, 
it will affect PLSs where the emission would otherwise orig- 
inate from the bulk of the shocked fluid (PLSs D, G, and 
E, which are discussed in §[2]3]). In this case, their appear- 
ance would resemble those of the fast cooling PLSs, since 
in both cases the emission arises from a very thin layer just 
behind the shock, and the image would become extremely 
limb brightened. This is potentially testable in a microlens- 
ing event or if a particularly nearby afterglow image will be 
well resolved. 

Finally, the afterglow image in general depends also on 
the magnetic field structure (its orientation) in the shocked 
region, not only on its absolute value (or strength). In this 
work it was assumed to be tangled on small scales with an 
isotropic distribution in the comoving frame of the emitting 
shocked fluid. In this case the emission and absorption co- 
efficients are also isotropic in that frame, which simplifies 
the calculation of the afterglow image, and it is very use- 
ful in this respect. However, such a magnetic field structure 
predicts no polarization of the afterglow emission, which is 
inconsistent with the measurement of linear polarization at 
the level of a few percent that has been detec ted in the op- 
tical or NIR afterglow of several GRBs (see ICovinol |2004| . 
and references therein). Therefore, in a future work we will 
discuss the changes that arise for other possible magnetic 
field structures in the emitting region. 
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3 DISCUSSION 

Analytic expressions were derived for the afterglow image for 
PLSs in which the emission originates from a very thin layer 
just behind the shock (§§[2II1[12J while simple semi-analytic 
expressions were obtained for the remaining PLSs in which 
the emission arises from the bulk of the shocked fiuid (ij |2.3|) . 
These expressions are for a rather general power law exter- 
nal density profile, pcxt oc r~'^ with fc < 4, for which the fiow 
is described by the BM76 self-similar solution. The relevant 
expressions are given in §[21 and illustrated in Figures [T] and 
[21 These expressions fully ag ree with the aft e rglow images 
that were shown and used in lGranot fc Loebl (|200ll ). which 
were calculated using the formalism of GS02. The fiux nor- 
malization for the different PLSs, which also provides the 
surface brightness normalization, can be found in Table 1 of 
GS02. 

The magnetic field in the shocked external medium 
must be considerably amplified at the shock in order to re- 
produce the afterglow observations. However, it is not yet 
clear how far downstream this shock produced magnetic field 
persists. It could in principle decay considerably at some fi- 
nite distance behind the shock. If the magnetic field is sig- 

This is obtained since for k{b -|- 1) < —(136 -I- 1) the integral 
in Eq. I I25I I becomes dominated by a narrow range of y values 
near the lower limit of integration, where y ~ j/min ~ C^x^ and 
X ~ (.y-yn-,in)/y'^in i° ^'^^ ''^'^g^ (CfcX^)^-* <y-y^in < CkX^, 
and the lower limit of integration dominates. 
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